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We compute the charm quark mass in lattice QCD and compare different formulations of the heavy quark, 
and quenched data to that with dynamical sea quarks. We take the continuum limit of the quenched data by 
extrapolating from three different lattice spacings, and compare to data with two flavours of dynamical sea quarks 
with a mass around the strange at the coarsest lattice spacing. Both the FNAL and ALPHA formalism are used. 
We find the different heavy quark formulations have the same continuum limit in the quenched approximation, 
and limited evidence that this approximation overestimates the charm quark mass. 



1. INTRODUCTION 

The charm quark mass is a relatively unknown 
quantity, yet this is a quantity which is read- 
ily calculable on the lattice. The Particle Data 
Group 1 quote 



1.0 < m 



charm 



(m, 



charm 



< 1.4 GeV 



(1) 



Most lattice calculations to date have been done 
in the quenched approximation. In this work 
we compare the results of the quenched contin- 
uum limit from three lattice spacings, using dif- 
ferent heavy quark formulations to results with 
two flavours of dynamical fermions with a mass 
around that of the strange quark. An idea of how 
heavy the sea quarks are can be seen by compar- 
ing the value of mps/my ~ 0.7 to the physical 
value of m T /m p ~ 0.18. Preliminary results were 
presented at the previous lattice conference [2] 
and results for the D s spectrum were presented 
in 0. 

2. DETAILS OF THE CALCULATION 

We used the Wilson plaquette gauge action and 
the Shcikholcslami-Wohlert action with a non- 
pcrturbative (NP) determination of the coeffi- 
cient csw- At each lattice spacing we computed 
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Table 1 

Ensembles of gauge configurations. K sea = de- 
notes a quenched ensemble. Nq denote the num- 
ber of configurations. 

09, Ksca) Volume a" 1 (GeV) N G 



(6.2,0) 24 3 x48 

(6.0, 0) 16 3 x 48 

(5.93,0) 16 3 x 32 

(5.2,0.1350) 16 3 x 32 



2.93 
2.12 
1.86 
1.88 



216 
302 
278 
395 



different definitions of the quark mass, for differ- 
ent values of the hopping parameter. We fixed 
the charm mass using the tuds pseudoscalar me- 
son. The details of the matching between the 
quenched and dynamical ensembles can be found 
in 0]. The details of the ensembles are shown in 
Tabled Throughout this work the scale is set by 
ro = 0.5 fm. 

2.1. Definitions of quark mass 

The quark mass can be defined in a number of 
ways. It can be defined from the quark coupling, 
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where K cr ;t is the value of the hopping parame- 
ter which corresponds to zero quark mass. The 
renormalised quark mass, in the mass indepen- 
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dent scheme (ALPHA), is then 

amy = Zm (1 + b m amo) amo (3) 

Another quark mass can be defined from the 
Axial Ward identity, and correlation functions 



2£^(z)^(0)) 
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where A is the improved Axial quark current, P 
the Pseudoscalar density, ans i and j label quark 
flavour. The subscripts i,j denote different quark 
flavours. For heavy-light correlation functions, 
the renormalised quark mass is 
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- (f + (b A - bp)arrio) am LL ] 



where mo = ^(itiq + TOq ). 

The values of all the improvement and renor- 
malisation coefficients have been determined to 
one loop in tadpole improved perturbation the- 
ory [H] in the MS scheme, as the NP determi- 
nations are not available for the coarsest lattice 
spacings. 

We also use the effective field theory formal- 
ism known as the FNAL method. This presumes 
the dominance of 0(am) n lattice artefacts over 
O(ap) 2 . This can be seen by considering the dis- 
tortion of the dispersion relation by the lattice, 

= Ml + §V + 0(/) (6) 

M 2 

where M\ is the rest mass and M 2 is the kinetic 
mass, defined by 

f _ d 2 E 
~M 2 ~~dpl 

It is the kinetic mass which is important for the 
dynamics of states containing a heavy quark. The 
definitions of quark mass have amo corrections to 
all orders in a, to a particular order in a$- At 
tree-level the quark mass is, 
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ami — log(f + amo) 



(8) 



A perturbative definition of the kinetic mass is 
give by 

iroi s i n h(ami) 



am 2 (ami) 



1 + sinh(ami) 



(9) 



These quark masses can be used to determine the 
kinetic hadron mass, 



aM 2 T — aM\ + (am 2 — ami) 



(10) 



The tree-level quark masses give a poor estimate 
of the kinetic mass shift. However, using the one- 
loop masses, the perturbative M 2 tracks the NP 
M 2 from the dispersion relation, without the in- 
crease in statistical errors coming from non-zero 
momentum states. The one-loop expressions re- 
lating the lattice quark mass and the pole quark 
mass are, 



m^ = rrii + 9 2 ^ml tanhm^ 
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where the Z functions are given in 0. We now 
refer to m[\ as mi, 2 for brevity. These pole quark 
masses are then matched to the MS mass using 
the one-loop expression from [S]. 

We interpolate each of the four quark masses 
{mji,my ,mi,m 2 } from equations (|3*|).ll^|l. (|ll|) 
and 11 2|) against the pseudoscalar heavy-light me- 
son mass, where the light quark mass has been 
interpolated to the strange quark mass. We de- 
note this Mhs ■ AH the quark masses have correc- 
tions which scale with the quark mass, but the 
inverse hadron mass scales with quark mass, so it 
is unclear whether to plot amQ vs. 1/Mhs or vs. 
Mhs- We do the latter as we are interpolating 
in a finite range of amQ, where a polynomial in 
amq can be expanded in terms of 1/amQ. The 
heavy quark interpolation for m A on the matched 
ensembles is shown in figure ^ For ensembles 
with four different heavy quarks (/3 = {6.2, 6.0}) 
we use a quadratic function, for ensembles with 
three (f3 — {5.93.5.2} we use a linear function. 

The resulting charm mass is in the MS scheme 
at the scale of the lattice spacing. We evolve this 
mass using the renormalisation group equation 
evaluated in a package called RunDec [H] to the 
scale invariant mass. We can then compare the 
different quark masses and take the continuum 
limit of the quenched data. 

2.2. Comparing quark masses 

Shown in figure [21 is the quenched continuum 
limit. The leading lattice artefacts are O(asam), 
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Figure 1. The hadron mass versus the itia quark 
mass for the ensembles (/3 = 5.2, n = 0.1350) and 
P = 5.93, k = 0). 
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Figure 2. The quenched continuum limit of the 
scale invariant charm quark mass. The closed 
(open) symbols show the quenched (dynamical 
fcrmion) data. The shaded circle is the contin- 
uum limit of the ALPHA masses and the diamond 
(offset for clarity) the FNAL masses. 



0(asa) and 0(a 2 ). For the masses itia and my 
these are clearly large and a linear extrapolation 
is not possible. However, they should have a com- 
mon continuum limit, as demonstrated by Rolf 
and Sint [o]. We enforce the common continuum 
limit and perform a simultaneous quadratic fit to 
both masses. Within large errors this is consistent 
with the result from Rolf and Sint. We perform a 
similar fit to the FNAL masses. The two methods 
agree in the continuum limit. 

The open symbols show the dynamical data. 
The dynamical ALPHA masses are lower than the 
quenched data. However, this is not repeated for 
the masses mi and m2, and the sea quark mass 
is still rather heavy. Moreover, the dynamical 
and quenched data could be effected by different 
lattice artefacts. 

Our analysis has larger statistical and system- 
atic uncertainties than the quenched result of Rolf 
and Sint, and other calculations. However, we 
have demonstrated that the charm mass in the 
FNAL and ALPHA schemes has the same con- 
tinuum limit. We have also found evidence which 
suggests the mass of the charm quark is too high 
in the quenched approximation at fixed lattice 
spacing, even though the mass of the sea quarks 
in the dynamical calculation is still rather heavy. 
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